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Subroutine Allows Easy Computation in Extended Precision Arithmetic 


The problem: 

To devise a program that allows relatively simple 
computation of very large numbers or very small frac- 
tions with extreme accuracy. Ordinary double preci- 
sion FORTRAN arithmetic is limited to numbers be- 
tween 10r 38 and 10+ 38 , and to 16 digits of accuracy. 


The solution: 

A subroutine called NPREC that can perform 
mathematical operations in extended precision float- 
ing point arithmetic. This subroutine handles num- 
bers that consist of 35 binary bits (1 word) for the 
exponent and 70 bits (2 words) for the fraction. A 
programmer can perform a variety of mathematical 
operations by writing standard FORTRAN state- 
ments within the NPREC routine instead of calling a 
separate routine to perform each operation. 


How it's done: 

NPREC can be used on any of the IBM 709/7090/ 
7094 computers. A statement, CALL NPREC, is 
used to enter the routine and subsequent statements 
are interpreted by the routine until the exit statement, 
CALL NPOUT, is encountered. Any statement that 
is interpreted in NPREC will take longer to execute 
than if it were executed outside of the routine, and 
thus all normal single or double precision statements 
should be executed outside of the routine. 

The NPREC routine incorporates subroutines in 
converting to and from extended precision numbers and 
also includes the following extended precision ele- 
mentary functions: square root, sine, cosine, arc 
tangent, natural log, and exponential. The only “built- 
in” library functions that may be used in NPREC are 


ABSF, MAX1F, MIN IF, AND SIGNF. There is 
limited use of FORTRAN exponential notation in 
NPREC; and extended precision numbers may be 
raised to a fixed point power of 2, 3, 4, 5, or 6 only. 
Any other type of exponential operation must be 
performed using the exponential and natural log rou- 
tines within NPREC. No more than 15 operations 
should appear in any single statement in NPREC. 

The NPREC routine will detect 20 errors. When one 
of these errors is detected, a comment will be written 
in 1 of 2 formats on output tape 6, after which a core 
dump is given and the job is finished. 

The time it takes to perform an arithmetic operation 
in NPREC varies with the values and relative values 
of the operands, and also depends on whether or not 
the operands are stored in COMMON. As an indica- 
tion of the time requirements, an extended precision 
division takes approximately 1.4 milliseconds on the 
7090 and a multiplication approximately 0.6 milli- 
seconds. The natural log and exponential functions 
require approximately 36 and 48 milliseconds respec- 
tively. Multiplying the 7090 time by 6 gives the 709 
time. 

NPREC and its associated subprograms require 
2053 storage locations. 

Note: 

Inquiries concerning this innovation may be di- 
rected to: 

Technology Utilization Officer 
Marshall Space Flight Center 
Huntsville, Alabama 35812 
Reference: B66- 10504 
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This document was prepared under the sponsorship of the National 
Aeronautics and Space Administration. Neither the United States 
Government nor any person acting on behalf of the United States 


Government assumes any liability resulting from the use of the 
information contained in this document, or warrants that such use 
will be free from privately owned rights. 


Patent status: 

No patent action is contemplated by NASA. 

Source: R. L. Berggren and J. C. Gysbers 
of North American Aviation, Inc. 

under contract to 
Marshall Space Flight Center 
(M-FS-1 136) 


Brief 66- 10504 


Category 01 


